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Nitrogen-vacancy (NV) centers are defects in diamonds, which, due to their electronic structure, 
have been extensively studied as magnetic field sensors. Such held detection applications usually 
employ the NV centers to detect held components aligned with the direction of the internally-dehned 
spin axis of the NV center. In this work we detect magnetic helds which are slightly misaligned with 
the NV center axis. In particular, we demonstrate that the NV center can measure the square of the 
angle between the magnetic held and the NV center axis with high sensitivity which diverges as the 
external held approaches a value pre-dehned by NV center’s internal parameters, in agreement with 
predictions. These results show that NV centers could be used as sensitive transducers for making 
quantum nondemolition (QND) measurements on systems such as nanomechanical oscillators. 


I. INTRODUCTION 

A variety of hybrid systems involving mechanical os¬ 
cillators coupled with spin systems have recently been 
studied and used for characterization of quantum prop¬ 
erties of their components 1 3 I With the advent of such 
systems, a class of measurement known as quantum non¬ 
demolition (QND) is re-emerging. Since the introduc¬ 
tion of QND measurements^®, this measurement class 
has found use in a wide variety of applications such as 
gravitational-wave detection® weak force detection using 
harmonic oscillators 7 } opticd®, etc. QND measurements 
allow repeated quantum measurements of certain quan¬ 
tities separated in time to give easily understood results. 
In a previous work, we proposed a hybrid system com¬ 
prising of a nanomechanical oscillator (NMO) coupled to 
a nitrogen vacancy (NV) center in presence of an external 
magnetic held with a QND interaction between the NV 
center and NMO. The QND measurement results from 
the quadratic coupling between the NMO and the NV 
center, which was predicted to have high sensitivity. 

NV centers are naturally occurring defects in diamond, 
which consist of a substitutional nitrogen atom and an 
adjacent carbon vacancy. These defects can also be 
artificially introduced in diamonds through irridiation 
and annealing 10 11 . The as-formed neutral NV center 
acquires an additional electron from elsewhere in the lat¬ 
tice, thus having six electrons 10 . This energetically fa¬ 
vorable charged state of the NV (NV - ) will be referred 
to as NV throughout this work. Owing to their ease 
of manipulation and readout, and long coherence times, 
they have recently being extensively characterized and 
used in various applications 11 . One such application of 
NV centers is as a sensor for steady and oscillating mag¬ 
netic fields 12 17 . These previous works have shown that 
the NV centers can be used to detect nanoTesla mag¬ 
netic fields with nanoscale spatial resolution at tempera¬ 
tures below 4 kP or at room temperatur^®, and can be 
used to map out the magnetic fields in a two-dimensonal 
aresP 6 . Most of these applications choose the NV center 
aligned in the direction of the external magnetic field— 
or they detect the component of the magnetic field that is 


aligned in the direction of the NV centeiE®. Applications 
that detect magnetic fields that are misaligned with the 
direction of the NV center do so by iteratively fitting pa¬ 
rameters and minimizing error functions 1 ®. Due to these 
limitations, the maximum magnetic field measu red by 
such techniques has been limited to <10 mT 14 I®. 

In this work, we experimentally demonstrate that an 
NV center in an external magnetic field can be used as a 
detector with quadratic coupling to the angle 6 between 
the magnetic field and the NV center axis. Instead of 
simply considering the NV center to be a sensor for the 
magnetic field in a single direction, in this work we exper¬ 
imentally confirm a fuller but simple quantum treatment 
of the coupling which predicts that an NV center can 
be used to measure 0 2 with sensitivity which diverges as 
the magnetic field approaches a critical value, which may 
allow the NV center to be used for QND measurements 
in nano mechanical systems. We briefly summarize the 
previously reported theoretical results^ below. 

II. THEORY 

The energy structure of the ground state of an NV cen¬ 
ter has been extensively studied®®} and is understood to 
be as shown in Fig. [l] The ground state ( 3 A) is a spin 
triplet, with |0), |1) and | —1) spin sublevels, defined in 
the S z basis (z being along the NV center axis), where 
state \i) has m s = i. In the absence of any external mag¬ 
netic field, |1) and | —1) are degenerate, with a zero-field 
splitting of A = 2.87 GHz above the |0) state. In presence 
of an external magnetic field, the degeneracy between the 
|1) and | —1) states is lifted by Zeeman splitting. 

The Hamiltoni an fo r the ground state of the NV center 
can be written ad 13 * 18 !: 

H nv = hASz + g/iB (s • B e xt ) (1) 

where A is the zero-field splitting, g ^ 2 is the electron 
g-factor in the NV center, /ib is the Bohr magneton, S 
is the electronic spin of the NV center and B ext is the 
external magnetic field. 
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FIG. 1: Energy structure of the ground state of the NV 
center. 


In our experiment, the applied magnetic field B ext may 
not be aligned with the NV spin vector S. This gives the 
term S • B ext in Eq. [l]an angular ( 6 ) dependence, where 
6 is the angle between B ext and S. For this work, we 
fix the direction of the NV center, and vary the exter¬ 
nal magnetic field. For the NMO-NV proposed hybrid 
systenP, the direction of magnetic field is fixed and the 
direction of the NV center moves due to the oscillation of 
the NMO. In either case, H^v takes the following form, 
where we assume B ext is in the y — z plane: 


H n v = HAS 2 Z + g^BBext (Sy sin 6 + S z cos 0) (2) 


The eigenvalues of H^y, A^(0), with 6 as a parameter, 
have the form 


Xi(0) — huii ± KiO 2 ± O(0 4 ) (3) 

where i m —1, 0, or +1 is a label which is equivalent to 
the NV center m s when 6 = 0. TioJi is simply the result of 
combining the zero field splitting with the static external 
field at 0 = 0: 


u±i = A (1 ± B ext /B z f s ) (4) 

^o=0 (5) 


where B z f s = hA/gys ~ 102.5 mT is the magnitude of 
the effective internal magnetic field which results in the 
zero field splitting A. 

The 0 2 coefficient in A i(0), /^, is: 


k±i = —hA 


kq = hA 


B e 


2{B ext is: B zfs ) 

B 2 ext 


B 2 

ext 


B lfs 


(6) 

( 7 ) 


Therefore, the 6 2 dependence of the energy gap between 
|0) and | — 1), and |0) and |±1) near 6 = 0 (kq— and fto±, 
respectively) are 


k,q- = —hA 
/^o+ = —hA 


(B ext )(3B ext -\-B z f s ) 

2 ( B %s- B2 e X t ) 


( B ext ){3B ext — Bzfs) 


( 8 ) 

(9) 


The denominators of k,q — and fto± suggest that the sen¬ 
sitivity of the spin transition frequencies to 0 2 around 
6 = 0 should increase much faster than linearly and even 
diverge as B ext approaches B z f s , the result of an avoided 
level crossing in the NV center. 


III. EXPERIMENTAL DATA 

To test the predictions, electronic spin resonance 
(ESR) was performed on diamond nanocrystals with 
NV center ensembles. The nanocrystals used for this 
purpose were ^100 nm and were prepared using high 
pressure high temperature (HPHT) techniques (Adamas 
Nanotechnologies: ND-NV-lOOnm (400NV)). The dia¬ 
mond nanocrystals were spin-coated on polished silicon 
with low density, so individual nanocrystals could be re¬ 
solved. Confocal laser scanning microscopy (CLSM) was 
employed to visualize and detect the photoluminescence 
of the NV centers. Green laser light (A = 532 nm) was 
used to excite the NV centers in these nanocrystals and 
their photoluminescence was detected using a single pho¬ 
ton counting module (SPCM) after filtering out wave¬ 
lengths between 517 nm and 548 nm using a notch filter. 
To drive resonant transitions between |0) and |d=l) spin 
states, a microwave (MW) drive was applied to this sam¬ 
ple using an un-terminated loop, 2 to 3 mm in diameter, 
made out of 30 fi m diameter gold-plated tungsten wire. 
A dip in the photoluminescence intensity is detectable 
when the NV centers are driven at a frequency resonant 
with a spin transition. In the absence of an external 
magnetic field a characteristic ESR signal, a dip in pho¬ 
toluminescence at ^2.87 GHz i s expect ed from single NV 
centers or NV center ensemble d 11 * 14 * 19 *, and was obtained 
on our samples as well. 

An external magnetic field B ext with continuously vari¬ 
able magnitude and direction was applied to the nanodi¬ 
amonds using three pairs of electromagnetic coils. Due 
to experimental limitations, ~ ±100 mT was available 
in the direction perpendicular to the sample (£'), and 
« ±25 mT in each of the two directions parallel to the 
sample surface ( x' and y'). The range of accessible mag¬ 
netic fields is illustrated graphically in Fig. 2(a). For 
weak fields, B ext < 25 mT, the external magnetic field 
could be applied in any direction, but for B ext approach¬ 
ing B z f s , only a small angular range was accessible. 

Individual diamond nanocrystals were randomly ori¬ 
ented on the sample. Within each nanocrystal, the di¬ 
rection of an NV center axis can be in any one of the 
four possible orientations along the carbon tetrahedron 
bonds, defining the NV center z direction (shown with 
blue arrows in Fig. 2(b)). Due to the symmetry between 
|1) and | —1) states, the NV z direction could also be cho¬ 
sen to be any of the four orientations directly opposite to 
the carbon tetrahedron bonds (shown in red arrows in in 
Fig. 2(b)). If each diamond nanocrystal has many NV 
centers with random orientations, there are eight possi¬ 
ble directions where B z f s and B ext could align with each 
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FIG. 2: Pictorial 3-D representation of experimentally accessible external magnetic field (B ext ) and possible NV 
center orientations (directions of B z f s ) in any arbitrary diamond nanocrystal (a) The pink sphere represents the 
direction of possible B z f s for all diamond nanocrystals. The cuboid with square faces of side 25 mT and a height of 
±100 mT represents the experimentally accessible magnetic field region, along the experimentally defined z' 
direction of B ext . This limited the angular range of a certain magnitude of magnetic field we could apply around the 
sample, as shown by the shapes inside the cuboid, (b) A carbon tetrahedron structure in any arbitrary diamond 
nanocrystal is shown. Assuming the nitrogen atom is at the center of the tetrahedron (pink), and one of the other 
carbon atom (green) sites is vacant (light green), the NV-axis could be directed along one of four possible 
orientations (shown in blue arrows). Due to the symmetry between |1) and | —1) spin states, the NV-axis could also 
be directed along the four orientations opposing the carbon tetrahedron bonds (shown in red arrows). 


other, each selecting a subset of the NV centers. These 
eight directions define the orientations along which we 
expect the NV centers could behave as 0 2 sensors with 
very high sensitivity. 

Experimentally, the first task was to identify the NV 
axis orientations in a particular diamond nanocrystal by 
finding the magnetic field directions at which the ex¬ 
pected Zeeman splitting is observed. While maintaining 
a fixed magnetic field magnitude, ESR spectra was mea¬ 
sured at discrete steps in magnetic field directions. The 
magnetic field magnitude was chosen to be weak enough 
that our coils could produce it in any direction, although 
due to the symmetry of the NV centers, sweeping the 
magnetic field only along one hemisphere was sufficient. 
At each field orientation, a microwave drive was swept in 
frequency, and dips in the florescence indicated resonance 
with a spin transition. 

We monitored the ESR dip depth and location while 
sweeping a fixed magnitude of B ext around in polar [0- 
9O°](0) and azimuthal [0-360°] (0) angles, where the po¬ 
lar axis is defined to be normal to the substrate. The 
results for a specific diamond nanocrystal, which was ~ 
10 fim away from the edge of the MW loop, are shown in 
Fig- m Similar behavior was observed in other diamond 
nanocrystals within 20-30 /i m of the MW loop. Fig. |3a| 
and Fig. [3b] show the behavior of the ESR dip depth and 
dip location respectively, while sweeping B ext in a hemi¬ 
spherical manner with 20 mT magnitude. Such plots help 
us broadly recognize the regions where the ESR dip depth 
is a maximum and dip location is a minimum, as shown 


by the circled regions in these plots. These regions give 
the approximate orientations of the NV centers in this di¬ 
amond nanocrystal, where we expect them to behave as 
0 2 sensors with very high sensitivity when the magnetic 
field magnitude is increased. Due to the experimentally 
accessible magnetic field region being close to the pole 
of the hemispheres, we chose the region nearest the pole 
for further exploration. Fig. [3d] shows the ESR dip in 
this region which fits well to a Lorentzian lineshape with 
frequency at ~ 2.3 GHz, as expected. 

We further narrowed down the spread of the NV center 
response in both 0 and 0 by probing with higher mag¬ 
netic fields. Fig. [3c] shows such a plot of the behavior of 
the ESR dip location while sweeping an 80 mT B ext . As 
can be seen in this plot, the spread of both polar and 
azimuthal angles narrows down as the external magnetic 
field is increased. In this case, the spread in the polar an¬ 
gle is much finer (<5°) than the spread in the azimuthal 
angle (<25°), because the field direction is near the pole. 
The shifted ESR dip at 80 mT at 13.4° </> and 267° 0 is 
shown in Fig. [3e] fitted to a Lorentzian lineshape with 
the dip frequency at ~ 0.63 GHz. 

By sweeping B ext between 0-95 mT in 0 and 0 space, 
the ESR dip was observed to shift down linearly in fre¬ 
quency on increasing the magnitude of B ext , as seen in 
Fig. [4a] Fitting this experimental data to a straight line, 
we get an ^-intercept of 102.46 mT, which is very close to 
the expected value of B z f s = 102.5 mT, and a ^/-intercept 
of 2.863 GHz, which is very close to zero-field splitting, 
A = 2.87 GHz. At several magnitudes of B ext , the vari- 
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FIG. 3: Locating NV orientations for an arbitrarily chosen diamond nanocrystal: © monitoring the ESR dip depth 
while sweeping 20 mT B ext around in azimuthal and polar angles about the +£ direction and projecting the 
hemisphere onto a polar plot and © similarly monitoring and plotting the ESR dip location. Circled regions in ([a]) 
and © exhibit maximum dip depth and minimum dip frequency, which can be thought of as orientations of the NV 
centers in a diamond nanocrystal. (|c| Similarly monitoring the ESR dip location at 80 mT B ext and projecting the 
spherical cap onto a planar plot. The colormap for the projection plots is also shown. Q ESR dip at 20 mT in a 
region closest to the pole, where the dip location reaches a minimum in frequency (12.5° (j) and 264° 0). 0 ESR 
dip at 80 mT at the same location as above, where the dip frequency reaches a minimum. 


ation of the ESR dip frequency, with respect to </> and 0, 
was found to be quadratic, as in Fig.|4b|and Fig. [4c] The 
coefficients of the 0 2 dependence at different magnitudes 
of B ext is shown as the sensitivity in Fig. [4d]and are in 
excellent agreement with the expected theoretical behav¬ 
ior. This experimentally verifies that the NV centers are 
highly sensitive 0 2 sensors, with sensitivity that increases 
as B ext increases, and diverges as it approaches B z f s . For 
comparison, if the NV center were only affected by the 
component of the magnetic field along its z axis, the 0 2 
sensitivity would be linear with B ext , as in Fig. [4d] 


IV. CONCLUSION 

In conclusion, we have shown that NV centers in a di¬ 
amond nanocrystal can be used as a high-sensitivity 0 2 
sensors. To use a particular nanocrystal for this purpose, 
a weak magnetic field is first swept over a wide range of 
angles to identify the orientations of the NV centers. Af¬ 
ter increasing the magnitude of the field and aligning it 
near the NV center z axis, the sensitivity of the NV cen¬ 
ter ESR frequency to 0 2 is greatly enhanced, even diverg¬ 
ing as B ext approaches B z f s ps 102.5 mT, as predicted in 
our previous work 3 . The ability of the NV centers to 
be highly sensitive 0 2 sensors is critical for their use in 
making QND measurements of e.g. the number state of 
a torsional harmonic oscillator. This class of measure¬ 
ments may be valuable in a variety of applications such 
as gravity wave detection and probing quantum behavior 
ofNMOs. 
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FIG. 4: Various trends seen with the shift of ESR dip on applying external magnetic field. As predicted ([a]) the ESR 
dip shifts down in frequency linearly with the B ext magnetic field and © the ESR dip depth and location have a 
quadratic behavior with angle; dip location is plotted here with respect to the polar angle at various magnetic fields. 

@ The predicted and measured absolute transition frequencies are in excellent agreement. © Plotting the 
quadratic sensitivity at different magnetic fields shows that the NV center is a much more sensitive 0 2 sensor than 
any L B Z sensor’, which measures the component of the magnetic field along its axis. 
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